ar X iv : p hy si cs / 0 51 21 30 v 1 1 4 D ec 2 00 5 Theory of solvation in polar nematics

We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The formal theory constructs the solvation response function from projections of the solvent dipolar susceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phase. Based on the properties of the solvent susceptibility from simulations and the formal solution, we have obtained a formula for the solvation free energy which incorporates experimentally available properties of nematics and the length of correlation between the dipoles in the liquid crystal. The theory provides a quantitative framework for analyzing the steady-state and time-resolved optical spectra and makes several experimentally testable predictions. The equilibrium free energy of solvation is anisotropic in the nematic phase and is given by a quadratic function of cosine of the angle between the solute dipole and the solvent nematic director. The sign of solvation anisotropy is determined by the sign of dielectric anisotropy of the solvent: solvation anisotropy is negative in solvents with positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous at the point of isotropic-nematic phase transition. The amplitude of this discontinuity is strongly affected by the size of the solute becoming less pronounced for larger solutes. The discontinuity itself and the magnitude of the splitting of the solvation free energy in the nematic phase are mostly affected by microscopic dipolar correlations in the nematic solvent. Illustrative calculations are presented for the Stokes shift and Stokes shift correlation function of coumarin-153 in 4-n-pentyl-4’-cyanobiphenyl (5CB) and 4,4-n-heptyl-cyanopiphenyl (7CB) solvents as a function of temperature in both the nematic and isotropic phase.